Analytical Results for the Statistical Distribution Related to Memoryless 
Deterministic Tourist Walk: Dimensionality Effect and Mean Field Models 
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Consider a medium characterized by N points whose coordinates are randomly generated by 
a uniform distribution along the edges of a unitary d-dimensional hypercube. A walker leaves 
from each point of this disordered medium and moves according to the deterministic rule to go 
to the nearest point which has not been visited in the preceding yi steps (deterministic tourist 
walk). Each trajectory generated by this dynamics has an initial non-periodic part of t steps 
(transient) and a final periodic part of p steps (attractor). The neighborhood rank probabilities 
are parameterized by the normalized incomplete beta function Id — J 1 / 4 [l/2, (d + l)/2]. The joint 

distribution Sj* d (t,p) is relevant, and the marginal distributions previously studied are particular 
cases. We show that, for the memory-less deterministic tourist walk in the euclidean space, this 
distribution is: S[™\t,p) = [T(l + I^)(t + I^ X )/Y(t + p + 1^)^,2, where t = 0, 1, 2, . . . , 00, T(z) 
is the gamma function and Sij is the Kronecker's delta. The mean field models are random link 
model, which corresponds to d — > 00, and random map model which, even for = 0, presents 
non-trivial cycle distribution [S£J? m (p) oc p" 1 ]: S^ m {t,p) = F(N)/ {F[N + 1 - {t + p)]N t+p }. The 
fundamental quantities are the number of explored points n e — t + p and Id- Although the obtained 
distributions are simple, they do not follow straightforwardly and they have been validated by 
numerical experiments. 



I. INTRODUCTION 

The random walks in regular and disordered media are 
a very explored subject, capable to model several phe- 
nomena, in particular the transport problems Q, Q • On 
the other hand, deterministic walks in disordered media 
are less explored and their behaviors are not completely 
understood yet. 

Here we consider deterministic walks in disordered me- 
dia. The disordered medium is characterized by N points 
whose coordinates are randomly generated, according to 
a uniform distribution, along the unitary edges of a d- 
dimensional hypercube. A walker leaves from each point 
of this disordered medium and moves according to the de- 
terministic rule to go to the nearest point which has not 
been visited in the last \i steps. The quantity \i is called 
memory and represents the required time (in number of 
steps) to the regeneration of the visited sites (refractory 
time). The trajectory generated by this dynamics (de- 
terministic tourist walk) has an initial non-periodic part 
of t steps, called transient, and ends in a stable cycle of 
period p steps (attractor) where the same sites are vis- 
ited in the same order. Although the dynamic may be 
simply stated, it has a complex behavior, with non-trivial 
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results for yt > 2 0, & H IS 0] This rule may be relaxed, 
allowing the walker to visit the nearest sites with greater 
probabilities than the furthest ones (stochastic tourist 
walk) 8, 9|. Other than suggesting a possible mechanism 
for migration |3( this walk has been applied to thesaurus 
characterization p| and monkey movimentation 0. 

Depending on the system dimensionality d and the 
memory /i, several situations may be considered. The 
simplest case is the lazy tourist (fi = 0). The walker 
remains caught in the initial site. Thus all sites of the 
system are attractors of unitary period. The joint dis- 
tribution is So,d{t)P) = ^t.o<5p,i- Although this situation 
is trivial, its extension to the stochastic tourist walk 
and for the random map model |l(J is interesting because 
of its analycity (glass transition in the former and non- 
trival cycle distribution in the latter). For /1 — 1, the 
tourist has to leave the site where he/she is and moves to 
the nearest site among the remainders. Here the tourist 
knows, at each step, only the nearest site from his/her 
present position, but he/she does not recall any of the 
sites visited. In this case, the trajectories always end in 
two sites, which are mutually nearest neighbors (couple). 
Therefore, the period distribution is Si t d(p) — <5p,2, but 
the transient distribution is not trivial and will be ad- 
dressed here. The stochastic tourist walk with /1 = 1 
has been investigated in Ref. 0. The cases [i > 2 (treat 
with the quantity r = /i — 1) drastically differ from the 
preceding ones. Even for fj, = 2, it is possible to obtain 
trajectories with long transients and periods. 

The higher the euclidean space dimensionality is, the 
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weaker are the correlations (as triangular inequality, for 
instance) among distances. When d — > oo, these corre- 
lations may be neglected and the distances between the 
sites may be considered as independent random variables. 
Only the back/forwards symmetry Dij — Dji is pre- 
served. This is the random link model (RLM) , originally 
proposed by Mezard and Parisi an( A lately explored 
by Percus and Martin |l2j. 

If the back/forwards distances are different, the sym- 
metry Di j = Dj i is not preserved and all the N(N — 1) 
distances are random independent variables. In this case, 
the only reminiscent characteristic related to previous 
model is the null distance from a site to itself (D i;i = 0, Vz 
and fi = 1). A variant of the asymmetric random link 
studied by Derrida and Flyvbjerg [ljj, called random 
map model (RMM), consists to eliminate the restriction 
Di^ — 0, allowing the distance from a site to itself not 
to vanish. This distance may represent a cost for the 
tourist remains in a certain site, for instance. Unlike 
previous models, even for /j, = 0, the dynamics for this 
model gives rise to a complex cycle period distribution. 
This walk corresponds to the mean field approximation 
for the Kauffman's networks [l3|. 

The main objective of this paper is to generalize the 
geometric distribution and analytically obtain the prob- 
ability joint distribution for transient time and attrac- 
tor period for a ji = 1 deterministic tourist walk. Also, 
we study how these distributions are affected by the di- 
mensionality and the border/finite-size effects. We stress 
that Id and the number of visited sites n e are the rel- 
evant quantities of the problem. To solve this problem 
we have adopted the following strategy. Firstly, the solu- 
tions have been obtained in the limiting dimensionalities, 
i.e., d — ► oo and d — 1, and then the distribution for fi- 
nite d has been inferred through a generalization of the 
geometrical distribution. All obtained results have been 
numerically validated. 

This article is divided as it follows: In Sec. [n] we 
present the parameterization of Cox's equation by Id- 
In Sec. IIIII we analytically determine the transient dis- 
tribution for the RLM for arbitrary N . An analogy 
to the geometric distribution is then established. The 
subsistence and capture probabilities through the walk 
are defined, which will be the standard interpretation to 
treat the considered models. We then focus on the one- 
dimensional (ID) systems. With a simple modification 
in the algebraic formulation, the transient distribution 
for the infinite medium is obtained. Finally, we gener- 
alize the obtained results to systems with arbitrary di- 
mensionalities and numerically show their validity. In 
Sec. IIVI the joint distribution of the transient time and 
attractors period for the RMM is obtained. Final con- 
siderations are presented in Sec. [V] where the role of Id 
and n e = t + p to the joint distributions is stressed. In 
Appendix lAl some special functions are recalled amd in 
Appendix [B] the cumulative distribution is obtained in 
terms of the subsistence probabilities. 



II. REFLEXIVE NEIGHBORS 

In a Poissonic process of dimensionality d, the proba- 
bility that an arbitrary event is the nrfi 1 nearest neighbor 

of its own rfi 1 nearest neighbor is given by Cox's equa- 
tion [3: 
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The medium dimensionality is implicitly considered in: 
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(denoted by p in Ref . |l4j ) which is written in terms of the 
normalized incomplete beta function (Eq. IA2I and IA3jl . 
As shown in Appendix the approximation is justified 
for d^$> 1. This is the relevant quantity that parameter- 
izes the transient time distributions studied here. This 
quantity suggests a characteristic dimensionality (do = 8) 
from which the effects related to the dimensionality of the 
system can be neglected. 

In particular, for m = n, Eq. ^ gives the n-order re- 
flexive neighbors probability (Dacey's equation 
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The 1-order reflexive neighbors probability represents the 
probability of finding the medium attractors for [i = 1 
walkers, and Eq. reduces to: 

1 



p (d) = 



1 



(4) 



and Pf> = 1/(2 - e~ d l % / 'y/nd/8), when d > 1. The 
quantity p[ d ^ can also be interpreted as the null tran- 
sient time probabilitySi.d(i = 0) = in the /i = 1 
tourist walk. The value of P-^ ranges from 3/4, for 
d = 0, and diminishes monotoly, converging exponen- 
tially to the asymptotic value 1/2, as d — > oo. 16] Notice 
that Eq. 0] produces results for any real value d, which 
opens the possibility to the interpretation to non-integer 
dimensionality (fractals). Curiously, the analytical con- 



tinuation for d = leads to P^ =3/4 and not P^ v> = 1 
as one would possibly expect. 
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III. TRANSIENT TIME DISTRIBUTION 

Here we obtain the treansient time distribution for 
fi = 1 tourist walk. We start obtaining analytical results 
for the RLM (since the site distances are independent 
random variables), then for ID systems, and finally the 
result is inferred to an arbitrary dimensionality systems. 
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A. Random Link Model 

The RLM represents a mean field model for the eu- 
clidean space high dimensionality limit. Consider JV 
independent continuous random variables X\, X2, ■ ■ ., 
Xjv, with probability density functions (pdf's) /i(;ri), 
f 2(2:2) ■> ■ ■ ., /n(xn), respectively. The random variable 

Y = min{A q , X 2l ■ ■ ■ , X^} pdf g(y ) can be determined 
as follow [jjj] : The condition for Y to assume a given 
value y is that at least one of the variables X\, X2, . . ., 
Xjq is equal to y and all of the remainders are greater or 

equal to y: g{y) = J^Zi Mv) H.f&i)=l dx jfj( x i)- 

If all variables X\, X2, . ■ ., Xn have the same pdf 
f(x), then: g(y) = Nf(y)[l - F(y)] N -\ where F(y) = 
dxf(x) is the f(x) cumulative function. In partic- 
ular, if f{x) is a uniform between and 1 then g(y) = 
JV (1 — y) , with y varying from to 1. 

In the high dimensionality limit, the distances are in- 
dependent random variables, but what is the pdf? In 
fact, it does not matter to the tourist walk. Let us see 
the reason. For a given random variable X, following 
the pdf f{x), one can define another random variable 

Y = h(X), with some aimed pdf g(Y) simply impos- 
ing the condition that their cumulative distributions are 
equal, G(Y) = F(X). If g(y) > 0, for all possible values 
of Y , then the cumulative function G is a bijection, and 
the definition Y = G~ 1 [F(X)] leads to h = G" 1 o F. 
Thus, if the distances D{ j are particular values of a ran- 
dom variable X, with pdf f(x), it is then possible to 
obtain the random variable Y — h(X) with uniform pdf 
g{y) — 1 in the interval [0; 1]. In this case, the cumulative 
function G is the identity function, so its inverse is the 
identity function too. Hence, Y — F(X) = dx f(x). 
As f(x) > 0,Vx, one has that Y is a strictly crescent 
function of X. Therefore, if the distances D\ and D2 
obey the relation D\ < D2 in the X metrics, their cor- 
respondent values in the metrics Y will also do. For the 
deterministic tourist walk, both metrics yield the same 
trajectories, since in each step, it is not important the 
distance to be run, but only the neighborhood rank. 

The transient time probability distribution for a deter- 
ministic walk, with fi = 1, can be obtained for the RLM 
by noting that: i) the distance matrix is symmetric, ii) 
the distance between the sites are independent random 
variables and uniformly distributed in the interval from 
to 1, and Hi) the walk distance decreases at each step, 
until the tourist enters a 2-cycle. 

The walk distance x\ in the first step is the minimum 
of the JV— 1 independent random variables. In the second 
step, only JV — 2 new independent distances are available 
(because of the metrics symmetry, the distance to the 
first site has already been evaluated). Thus, for a tra- 
jectory with a transient size £ < — 2, it is necessary 
and sufficient that Xt+i < x t < . . . < x 2 < x\ and that 
x t +2 > £t+i, where Xj represents the minimum of the dis- 
til 

tances to the JV — j sites not visited up to the j step. 
This leads to the following transient time distribution: 
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(•5) 



x t +l 



The Fig. [3] schema shows the calculation of the t + 2 
chained integrals of the Eq. [3] for some particular val- 
ues of t. Summing these terms (and multiplying the 
sum by the produtory of the Eq. [SJ) , one obtains generi- 



cally: s[ N r l (t) = a x l\)=2 a jlY?k=i a k I Ei=o a ^ where 
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[N — (i + 3)/2](i + 2) and E 
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(t + 3 - j)/2](t + 2 - j). Calling k = t + 2-j leads to: 
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N - (t + 3)/2 j^N - (k + l)/2 



(6) 



where t = 0, 1, 2, . . . , JV" — 2 and the transient time distri- 
bution of the RLM in the thermodynamic limit (JV S> 1) 
is 



t+l 

(t + 2)! 



(7) 



which leads to t = e — 2 and of = e(l — i). 

Fig. [3 shows the transient time distribution for some 
values of JV and the convergence to the thermodynamic 
limit. Observe that, in particular, (0) = 1/2 = p^ 00 -* 

(Eq. EJ. An immediate property is that: Si^j(l) = 1/3. 

To measure the finite size effect for a given JV, the mean 
square error: £ 1)W (JV) = ZZol^M ~ ^(t)] 2 /^ 
has been considered. Numerically, one verifies that, 
for JV > 10, E lfH is a power law 0: E ljrl (N) w 
0,08611A- 3 . 

The probability the tourist enters an attractor is not 
constant along the walk, but it can be analytically ex- 
pressed as a function of the trajectory step. Let us start 
with the thermodynamic limit. The cumulative distribu- 
tion is F^(t) = 1 - 1/ (t + 2)! and Eq.Qcan be alterna- 
tively written as: S$ (i) = [1 - F^ (t - 1)] (t + l)/(t + 
2) = (t+l)/[(t+l)!(< + 2)], where l/(J+l)! represents the 
probability of the walker subsisting (not entering an at- 
tractor) in the t initial steps, and (t+l) / (t+2) represents 
the probability of the walker being captured (entering an 

attractor) in the (t + l) th step. This leads to: S[°$(0) = 
1/2, S-$(l) = 1/2 x 2/3, 5^/(2) = 1/2 x 1/3 x 3/4, 
5^(3) = 1/2 x 1/3 x 1/4 x 4/5 • • • . Making an analogy 
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FIG. 1: Chained integral calculation of In each integra- 

tion step, the number of integrals is doubled. The summations 
are indicated in the bottom line. 
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FIG. 2: Transient time distribution for the random link 
model. From bottom to top, the curves refer respectively 
to N = 6, 12, 25, 50, 100, 200 and to limit N > 1. 



form S[ N r j(t) = [l_(JV-t-2)/{[JV-(t+l)/2-l](i + 

2 )}] IlLi(^ - k - !)/[(# - fc / 2 -!)(*+ 1)], which al- 
lows us to obtain the finite size correction factor for the 

subsistence probability: 



(10) 



In terms of the subsistence probabilities, Eq. can be 
rewritten as: S$(t) = [l-«g2(t + l)] n*=i C?) an d 

the cumulative distribution is F^(t) = 1 — ^=1 gi'rlO') 
can be determined according to the Appendix FBI Ob- 
serve that using F^(t), one can easily verify that the 
transient distribution for finite size is also normalized 



B. One-Dimensional Systems 



to the geometric distribution, Eq. [7| can also be written 
as: 

siy® = [i - «s } (?■)(« + 1)] ri > (8) 

i=i 

where the success and, consequently, the failure proba- 
bilities depend on the extraction stage. The subsistence 
th 

probability up to the j step is: 

A similar reasoning can be applied to finite size systems. 
With some manipulations, Eq.|S|can be expressed in the 



As the RLM represents the upper limit for the dimen- 
sionality (d — > oo ) in the euclidean space, the ID model 
represents the lower limit, which is the easiest to analyze 
for finite dimensionality systems. In this Section, the 
transient time distribution for ID systems in the thermo- 
dynamic limit is analytically obtained. The demonstra- 
tion considers a semi-infinite medium, which establishes 
a surprising connection with the RLM. With a simple 
modification in the previous model, the transient time 
distribution for the infinite medium can be obtained. 



1. Semi-infinite medium 

The analysis of the semi-infinite disordered medium 
aims to: i) develop the calculation kernel for the tran- 
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sient distribution in the infinite medium, ii) reveal a non- 
trivial equivalence between ID systems and the RLM, 
and Hi) evaluate the edge effect. 

The semi-infinite medium can be thought as a set of 
uncountable points randomly and uniformly distributed 
on a line-segment, with a mean density of r points per 
unitary length. Consider site sj placed at the origin of 
the line-segment. 

This model is described by a ID Poissonic process. 
The sites can be viewed as events, which occur randomly 
as time flows. The distances between consecutive sites 
follow an exponential pdf: f(x) = re~ rx for x > and 
fix) =0 otherwise. 

The transient time distribution for a deterministic walk 
of a tourist who leaves from the origin of the system, with 
[1 = 1, can be determined by noting that the run distance 
decreases each step. For a trajectory with a t transient 
steps, it is necessary and sufficient that Xt+i < Xt < 
. . . < x<x < x\ and Xt+2 > Xt+i, where Xj represents the 
run distance in the j th step. This leads to: 



S(t) 



poo H~l pxj — i 

/ d Xl re- rxl Y[ dx 

J0 -_ 2 Jo 



re 



Axt+ire rXt+2 



X t +1 



t + 1 

(t + 2)! 



(11) 



where the schema of the Fig. illustrates the calculation 
of theses integrals. 

This expression is equal to Eq. Q the thermodynamic 
limit for the RLM. It is also interesting to observe that 
the medium mean density (r points by an unitary length) 
does not interfere in the transient time distribution, it 
only affects the distance the tourist runs each step. 

It may seem strange, at first look; that two models that 
represent opposite limits, concerning the dimensionality, 
present the same expression for the transient time dis- 
tribution. This equivalence can be understood noting 
that although there exist strong correlations in the dis- 
tance matrix for ID systems, the distances that really 
matter for fi = 1 tourist are those ones between consec- 
utives sites, and in fact, these distances are independent 
variables. Numerical simulations have revealed that this 
equivalence does not hold for /i = 2, because in this case 
second neighbor distances are also important [l6| . 



2. Infinite Medium 

With a simple modification in the previous model, 
one is able to obtain analytically the transient time dis- 
tribution for a "left-and-right hand unlimited" medium. 
Consider the tourist leaves from site s\ of the infinite 
medium. For a transient time t, as mentioned before, 
the only additional condition is that the distance xo 
(between the sites sq and s±) must be greater than the 
distance x\ (between the sites s\ and s-f). One must 
also multiply the expression by 2, since the walk can 
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3(2) 



2-3-4-5 



4/120 
5(3) 



FIG. 3: Calculation of chained integrals of Eq. 1111 Each 
integration step doubles the number of integrals. 



start either to the left or righthand side. Applying 
these cosiderations to the Eq. ^2 one has: S[™\t) = 
2 J °° daqre"" 11 /~ dx r e - rxo U'jLl fa*'* dxjre^™* x 
f Xt t dx t +2 r e~ rXt+2 ■ The integral in x is not linked 
with the other ones, and can be simply calculated. 
In this way, the schema of the Fig. [3] became a little 
different. In the last integration stage, the term — rx 
is added to the e~ rx , e~ 2rx , e~ 3rx . . . exponents 
Consequently, the denominators factors 2, 3, 4, ... 
from the last integral will be increased by a unity 
Hence: S[^\o) = 2x2/6, ^(1) = 2 x 3/24 
S[™\2) = 2 x 4/120, fig (3) = 2 x 5/720 

S'^ ) (4) = 2 x 6/5040, S$~ } (5) = 2 x 7/40320, .. 
. Generically, for t = 0, 1, 2, ... , oo, one has: 



2(t + 2) 
(t + 3)! 



(12) 
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In this way, we have analytically obtained the normalized 
transient time distribution for ID [i = 1 systems with 
N > 1. This leads to: t = 2e - 5 and of = 2e(2 - t) - 8 
. Notice that S[^(0) = 2/3 = Pf 5 (Eq-gJ). 

Similarly to the RLM, the probability the tourist is 
trapped in an attractor at each walk step in the ID case 
can also be analytically obtained. The transient time 
cumulative distribution is: Pfi (i) = Y?k=Q S[°^ (k) — 
1-2/0+3)! and the recursive form of Eq.dis s[™\t) = 
[1 - (t - 1)] (t + 2)/{t + 3) = 2(t + 2)/[(t + 2)!(t + 3)] . 
So that: S[™\0) = 2/3,5^(1) = 1/3 x 3/4,5^(2) = 

1/3 x 1/4 x 4/5, Sffifi) = 1/3 x 1/4 x 1/5 x 5/6. In this 
way: 

S[f{t) = [l-q%\t + l)]f[q<$>V), (13) 

J'=l 

where the success and failure probabilities depend on the 
stage of the extraction: 



«i,iU) = 



1 



(14) 



represents the subsistence probability (not getting an at- 
tractor) up to the j th step. 

The notable regularity in RLM and ID models makes 
us to consider the subsistence probability as the proper 
quantity to the generalization for arbitrary dimensional- 
ity systems. Contrasting to the RLM, we have numer- 
ically verified, for ID systems, that the transient time 
distribution has a weak dependence on N. [Tfil | 



C. Arbitrary Dimensionality 

Here we present arguments which allow us to predict 
the analytical form of the transient time distribution for 
the \i = 1 tourist walk in systems with an arbitrary di- 
mensionality for N ^> 1. As shown in Fig. 21 numerical 
simulation results have revealed that the transient time 
distributions for arbitrary dimensionalities lies between 
the analytically obtained limiting distributions (d = 1 
and d — > oo). 

Comparing Eas. 1131 and 151 one notices the same math- 
ematical structure. Further, the null transient trajec- 
tory probability in a d-dimensional medium is (Eq. [3J| : 



S[™\o) = P x w = 1/(1 + 7 d ). This is also the probability 
the tourist is captured in the first step. Thus, the subsis- 



>(<*) 



tence probability in the first step is g^^(l) = 1 
1/(1 + / £ / 1 ). Comparing Eas.l51and ll4l we infer that the 
subsistence probability for each step of the trajectory for 
arbitrary dimensionalities is given by: 




FIG. 4: Effect of the dimensionality on the transient distri- 
bution. From bottom up, the curves refer to d — 1, 2, 3, 5 and 



From these subsistence probabilities, it is possible to 
build a closed analytical expression for the transient 
time distribution for arbitrary dimensionality through 
the analogy with the geometric distribution: 



(*) 



r(t + 2 + p 



(16) 



3 + I d 



(15) 



with t — 0, 1, 2, ... , oo, which leads to t = 
e [r(l + p - T(l + p 1)] . Although the obtained 
results have been based on a conjecture, numerical sim- 
ulations have confirmed their validity. 

Consider the following remarks. One can write this 
distribution in terms of Id and in terms of the number 
n e = t+p of explored sites, once p = 2 remaining a t. The 
cumulative distribution can be obtained according to the 
AppendixEl F^\t) = 1 - T(l + p/T(t + 2 + p. 



IV. JOINT DISTRIBUTION IN THE RANDOM 
MAP MODEL 

The Derrida-Flyvbjerg random map |lOj is a mean field 
approximation for the Kauffman |l3j network. The map 
is built associating to each one of the N sites a random 
site; and the movement rule is to go, at each time step, 
to the successor site. Eventually a site may be its own 
successor. Even for fj, = and \x = 1 the model presents 
a non-trivial period distribution. This model may be 
applied to situations where the concept of distance rep- 
resents a cost, and the links are represented by a directed 
graph. In the RMM, contrasting to the preceding studied 
cases, even with [i = 0, it is possible to obtain periods 
varying from 1 to N. Therefore, a joint distribution for 
the transient time and attractors period is required to 
completely describe the walk. 
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Consider that the tourist starts the walk at site si. 
For the walk to have a transient t = and a period 
p = 1 (i.e., consists on a single point), the site fol- 
lowing s\ must be s\ itself. Hence, the probability is 

•So^M) = 1 / N - For the walk to have t = 1 and 
p = 1, the tourist must go to any site S2 among the 
N — 1 reminders, and in the following step to remain at 
a 2 , leading to: sg^l, 1) — (N — 1)/N x 1/N. Thus, 
the transient time marginal distribution for p = 1 is 
S$l(t, 1) = (N-l)/Nx(N-2)/N-. ■ (N-t)/Nxl/N. 
To obtain a t = transient and a p = 2 period, 
the walker must go to any of N — 1 remaining sites 
and, in the second step, return to s\\ Sg^^O, 2) = 
(TV - 1)/N x 1/N. For t = 1 and p = 2, the walker 
must go to any N — 1 remaining site, in the second step 
go to any of the N — 2 reminders, and finally return to S2'. 
5^(1,2) = (TV-l)/TVx(TV-2)/TVxl/TV. The transient 
time marginal distribution for p = 2 is: £q^(£, 1) = 
(TV- 1)/N x (TV-2)/TV ■ ■ ■ (N-t)/N x (N-t- 1)/N x 1/N. 
Generalizing this procedure for an arbitrary period p, one 
obtains the values of So^fap) displayed in the Tabled 
This A^-order matrix is symmetric and all elements below 
secundary diagonal are null. 



p 


t = 


t = 1 


t = 2 


1 


1 

w 


iV — 1 l 


W-l TV— 2 _^ 


2 


N — 1 _i_ 


W-l /V— 2 1 


TV-1 TV- 2 TV- 3 










3 


N — l TV — 2 l 

TV TV TV 


TV-1 JV-2 7V-3 1 
N N N N 


TV-1 TV-2 TV-3 TV-4 1 
TV TV TV TV TV 



TABLE I: Joint distribution for transient time and cycle pe- 
riods. The table is symmetric and the relevant quantity is the 
number of explored sites n e — t + p. 

With a simple inspection of the values in Table [I] one 
concludes that the transient time and attractor period 
joint distribution is: 



1 -pr N-j 



N 
T(N) 



r(TV-t-p+l)TV*+P 



(17) 



p and the transient t are identical, when time is retarded 
by one unity. Also, observe the strong influence of TV 
in the form of the distribution, which diverges in the 
thermodynamics limit. 

In terms of the number of explored sites n e = t + p = 
, N, Eg. 1 171 becomes: 



1,2, 



°Q,rm\ n e) 



n e T(N) 



t(tv- 



1)TV 



(18) 



Fig. El shows the validation of Eq. El through numerical 
simulation. 



2,0x10 



1,5x10" 



CO 



1,0x10' 



5,0x10 




500 



n 



FIG. 5: Distribution of the number of explored sites given by 
Eq. ll8l and numerical simulation. From the top to the bottom, 
the curves refers respectively to N = 1000, 2000, 4000, 8000, 
16000, 32000, 64000 and 128000 points by map. The greater 
N is, the wider the distribution becomes. 



The analogy we have estabilished to the geometric dis- 
tribution may also be applied to the RMM. Adopt as a 
failure the exploration of a new site and as success the re- 
visit to a given site. Consider that the walker leaves from 
site si. For the walker to explore the map in the first time 
step (visit a new site), site s\ must be connected to some 
of the TV— 1 other sites, but not to itself. The exploration 
probability in the first step is q^il) = (TV— 1)/TV. The 
walker then goes to the site S2- To explore a new site in 
the second step, site S2 must be connected to some of 
the TV — 2 reminder sites, but neither to itself nor to si. 
Therefore q ( N r ] m (2) = (TV - 2)/TV. This reasoning may be 



where t = 0, 1, 2, . . . , TV - p. Using this expression, generalized for an arbitrary time step j < TV: 



one can obtain the marginal probability for the attrac- 



Thus: Stfllnip) = Er=o P5 o!™(i>P) 



N-p 



.rrn V 



tor period p 

^L p r(TV)/T(TV - j + l)N j . This result agrees the 
ones obtained in Refs. [1113, with TV > 1: S^Jp) = 
VvW/^dy e-y 2 / 2 = 7V?W)erfc(p/v / 2TV) where 
(2/3- tt/8)TV. 

^ (P) = 



1o,rm(j) — 1 



(19) 



p = ^nN/8 and a. 

Notice that Table [I] symmetry implies to sjf^. 



Thus, the Eg. 1181 is rewritten in the form: 

= I 1 - Cl>e)] II iSlij) ■ ( 2 °) 



5oroi(^ — 1)' ^' e '' ^ ne mar S ma l distributions for the period 



The equivalence between the Eqs. ^| and |20 is imme- 
diatly verified. Notice that the exploration probability 
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decreases in arithmetic progression each step of the tra- 
jectory, contrasting to the preceding models, where the 
subsistence probability decreases in harmonic progres- 
sion. Another notable difference is that, in the thermo- 
dynamic limit, the exploration probability is unitary. In 
this way, the trajectory may have an infinite transient 
time, what characterize the chaos. 

The cumulative distribution for the number of explored 
site n e may be obtained (according to the Appendix lB|l : 
4 N l(ne) = 1 - T(N)/T(N - n e )N n °. 

V. CONCLUSION 

We have obtained the joint distribution for the fi = 1 
tourist walk for an arbitrary dimension in an euclidean 
space and also in two mean field models. These distribu- 
tions are parameterized by I e i and the number of explored 
sites n e =t+p. The former has been first introduced by 
Cox in the context of spatial statistics. Except for the 
RLM, we have not succeeded in obtaining a closed ana- 
lytical form in the finite size regime for the joint distri- 
butions in an arbitrary dimension. Except for the RMM, 
these distributions show a fast (factorial) convergence to 
the attractors. It is an open question whether this be- 
havior remains or not valid when the tourist has short 
range memory. 
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The normalized incomplete beta function [l8j is de- 
fined as: 

^ a < 6 ) = C Ate-\i-t) b ~ l (A2) 

B{a.b) J a 

with Re(a) > e Re(6) > and the beta func- 
tion ^3 is defined as being the normalization factor 
of I x (a,b): B(a,b) = 5(6, a) = dt ^(l - t) b ~ l = 
r(a)r(6)/r(a + b). One can conclude that the beta func- 
tion is an extension of the inverse of the combination of 
the Newton's binomial. 

It is also convenient to define the complementary func- 
tions: The complementar non-normalized gamma func- 
tion is written as: T(a, b) — f£° dt t a ~ 1 e~ t = T(a) — 
7 (a, b) and the complementar error function 

El 

is de- 
fined by: erfc(z) = 2 J°° dt e~* / 'y/n = 1 — crf(z). For 
\z\ 3> 1, the complementar error function has the follow- 
ing assymptotic form: erfc(z) w e~ 2 /(zy/n). 

Consider the behavior of Eq. IA2I when b 3> a ~ 
1. In this case, Eq. IA2I is written as: I z (a,b) « 
b a /* dt t°- x (l - t) b /T{a). If t < 1, then the factor (1 - 
tf = e 61 ^ 1 "*) w e - fct , so that: I z (a,b) w j(a,bz)/T(a), 
where 7(0,6) is given by Eq. IA1I When a = 1/2 
and b > 1 one has: I z (a,b) « j(l/2,bz)/T(l/2) » 
erf(\/&i) = 1 — cric(Vbz) and 

I,(o, 6) «1- -L e - 1 L. (A3) 
V w Vbz 
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APPENDIX A: SOME SPECIAL FUNCTIONS 

Some special functions which are extensively used 
along the text are reviewed here. Initially, consider 
the gamma function 18] T(z) = f™dt t z - x e~\ which 
main property T(z) = (z — l)T(z — 1) enables it to be 
a generalization of the factorial T(z) = (z — 1)!. For 
\z\ » 1 e |arg z\ < ir, Stirling's approximation ^3 is 
usedT(z) « v / 2Tr/z(z/e) z [l + l/(12z) + l/(288z 2 ) + •••]. 
The non-normalized incomplete gamma function |T3 | is 
defined as: 

,6 

7(0,6) = / det a_1 e- t (Al) 
Jo 

and presents the following property E3 7(1/2,2:) — 
2 Jq dt e - * 2 = \fir erf(y/x), where the error func- 
tion E3 1S defined as: erf(z) = 2/y / 7rJ 2 di e - * 2 = 

2 /V^J2kLo(- 1 ) kz2k+1 /[ kl ( 2k + !)]> which is monotoly 
increasing from erf(0) = to erf (00) = 1. 



APPENDIX B: CUMULATIVE DISTRIBUTION 
IN SUBSISTENCE PROBABILITY TERMS 

Here is convinient to deal with the cumulative of the 
generalized geometric distribution: 

t 

S(t) = [l-q(t + l)]l[q(j). 

3=1 

It can be obtained as follow: 

5(0) = 1-5(1) 

5(1) = g(l) - g(l)g(2) 

S(2) = g(l)g(2)-g(l)g(2)g(3) 



S(t) = q(l)q(2)-..q(t)-q(l)q(2)-.-q(t + l). 

Summing both members, yields: 

t t+i 
F(t)=Y,S(k) = l-l[q(j). (Bl) 

fe=0 j=l 
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